Speed of ion trap quantum information processors 



o : 
o . 
o ■ 

(N ; 

5h ! 
Oh- 

<: 
in : 

(N . 
(N ' 

> : 

. 

oo ■ 
o ■ 
cn : 
o . 
o ■ 

o ■ 
^ : 

^— > ■ 
c ■ 

3 : 



A. Steane*, C. F. Roos, D. Stevens*, A. Mundt, D. Leibfried, F. Schmidt-Kaler, R. Blatt 
*Oxford Centre for Quantum Computation, Department of Physics, University of Oxford, 
Clarendon Laboratory, Parks Road, Oxford 0X1 3PU, England. 
Institut fur Experimentalphysik, Universitat Innsbruck, Austria 
(February 1, 2008) 

We investigate theoretically the speed limit of quantum gate operations for ion trap quantum 
information processors. The proposed methods use laser pulses for quantum gates which entangle the 
electronic and vibrational degrees of freedom of the trapped ions. Two of these methods are studied 
in detail and for both of them the speed is limited by a combination of the recoil frequency of the 
relevant electronic transition, and the vibrational frequency in the trap. We have experimentally 
studied the gate operations below and above this speed limit. In the latter case, the fidelity is 
reduced, in agreement with our theoretical findings. 



I. INTRODUCTION 

Experimental methods which allow both coherent con- 
trol, and rapid and reliable measurement of the quantum 
state, have been available for some time for single ions 
held in ion traps. In an influential paper, Cirac and Zoller 
jlj showed furthermore how laser manipulation of a string 
of ions in a linear ion trap can achieve general coherent 
evolution of the joint state of the ions, using currently 
available technologies, and with good scaling properties. 
Specifically, the coherent evolution can be driven so as 
to realize any unitary transformation on the joint state 
of the ions, including transitions from product states to 
entangled states. Soon after their theoretical work, the 
essential ingredients of the method were realized experi- 
mentally for a single trapped ion Q), and more recently 
two ions were driven from a product state to an entangled 
state with around 90% fidelity using closely related ideas 
H . This is the only efficient creation of entanglement yet 
observed in any area of physics. 

The combination of universal driven unitary evolution, 
an exponential scaling of the available Hilbert space with 
system size, and reliable measurement of the resulting 
state, are the essential ingredients for a future quantum 
computer. The scaling of system size in this argument 
is measured by the way the whole physical apparatus 
becomes larger, and operation of quantum logic gates 
slower, as more qubits are added to the system. In the 
case of the linear ion trap, the adding of further qubits 
(ions) into the trap is straightforward. The growth of the 
physical apparatus, and the slowing down of the process- 
ing, is dominated by the optical and electronic equipment 
needed to control the ions. This growth and slow-down 
is a polynomial, not exponential, function of the number 
of qubits, for small numbers of qubits. Technical prob- 
lems will place a limit, as yet unknown, on the highest 
number of qubits which might be feasible, but below this 
limit the ion trap has the good scaling properties which 
allow the essential principles of quantum computing to 
be experimentally realised. Hence, laser-cooled ion trap 
experiments have a dual interest both for studying new 



avenues in atomic physics (e.g. interferometry with en- 
tangled particles) and for understanding quantum infor- 
mation processing in a specific realisable system. 

In this work we examine, both theoretically and exper- 
imentally, the speed with which general processing oper- 
ations can be driven in the ion trap system. We are 
concerned with the intrinsic limitations imposed by the 
physics of the system, such as the tightness of the trap 
and the presence of a rich energy level structure in the vi- 
brational modes. The initial proposal of Cirac and Zoller 
showed how to drive quantum operations in the limit of 
small Rabi frequency (slow operations), and initial exper- 
iments have been in this limit. Our aim here is to clarify 
the trade-off between precision and speed of the quan- 
tum logic gates, and to derive the optimal way to operate 
the processor. This optimum is a compromise between 
speed-associated problems such as off-resonant excitation 
of unwanted transitions, and the basic decoherence rate 
due to environmental coupling which will limit the per- 
formance when the gates are too slow. Our main observa- 
tion is to confirm the statement in Q| , that the maximum 
gate rate obtained by the Cirac Zoller method scales as 
the geometric mean of the trap vibrational frequency and 
the recoil frequency associated with the multi-ion string. 
This is in contrast to statements made elsewhere that the 
maximum gate rate is roughly proportional to the vibra- 
tional frequency, though this faster rate may be available 
if further tricks are adopted, such as excitation in the 
node of a laser standing wave Q . We extend the discus- 
sion in Q by making a quantitative statement, equation 
(pp|), of the relation between speed and precision of the 
gate. 

We also calculate the maximum gate rate for the 
method proposed by Monroe et al. Q , where a controlled- 
not operation is achieved without the need for an addi- 
tional transition and laser frequency, by driving multiple 
Rabi cycles on the carrier, for well-chosen values of the 
Lamb-Dicke parameter. 

We experimentally study the Cirac-Zoller method, us- 
ing a single trapped Calcium ion cooled to the ground 
state of the motion in one dimension. The experiments 
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are akin to those in [[2|]6|] where the methods were first 
demonstrated for one ion, except that we examine the 
regime as yet unexplored, where the gate rate is above 
the recoil frequency, so that the calculations carried out 
in the small Rabi frequency limit are no longer appropri- 
ate. 

The paper is set out as follows. Section || presents 
the theory of the logic gates, and calculations of their 
fidelity as a function of the main parameters, primarily 
the Rabi frequency of the atom-light interaction, and the 
trap tightness. The calculations involve numerical solu- 
tion of the Schrodinger equation for the system. We note 
the need to tune the laser light to the resonant frequency 
correctly, that is, taking into account the light shift (a.c. 
Stark shift) of the levels. We present results in the first 
instance considering only one mode of vibration of the 
ion string, and then in section II G we briefly discuss the 
influence of the other modes. Section III presents our ex- 
periments. We cool a single trapped calcium ion to the 
ground state of the motion in one dimension, and then 
observe Rabi flopping on the first blue motional sideband 
of the narrow 729 nm transition. By driving the sideband 
in the regime where the Rabi frequency is of the order of 
the trap vibrational frequency, we observe the expected 
trade-off between speed and precision of the operations, 
caused by off-resonant excitation of unwanted motional 
states. 



II. THEORY OF SWITCHING RATE 

A. Preliminaries 

Consider a string of two or more ions in a linear trap. 
The trap is strongly confining along x and y directions, 



and less strongly confining along the z axis. We will 
assume motion in x and y directions is unexcited, and 
consider one normal mode of vibration of the ion string 
along z, writing its angular frequency lu z . The other 
normal modes along z will be assumed to be unexcited 
throughout. (We will consider the element of approxi- 
mation introduced by this assumption at the end.) For 
brevity, we will refer to this vibrational degree of free- 
dom as "the normal mode" or "the vibration" . It has an 
evenly spaced ladder of energy levels E n = (n+l/2)hu z . 
Typically in an experiment one would choose the mode 
of interest to be the second or third, i.e. lu z = y3ui z cm 



y/W/Zl 



where i 



is the frequency of the cen- 



tre of mass mode. For modes other than the centre of 
mass, the Lamb-Dickc parameter describing the ion-light 
coupling will vary from one ion to another |Q, but this 
presents no problem as long as it is taken into account 
when choosing laser pulse intensities and/or durations. 
We will assume that the laser light is directed onto one 
ion at a time. 

The Lamb-Dicke parameter is rj = rji/^/N where N 
is the number of ions, and 771 = yj En/hu> z . Er = 
(rhk z ) 2 / (2M) is the recoil energy for a single ion ini- 
tially at rest undergoing a tt pulse interaction with the 
laser field; k z is the z component of the wavevector k of 
the laser field, M is the mass of the ion, and r = 1 for 
a single-photon transition, r = 2 for a Raman transition 
assuming the geometry ki z — —k,2 Z for the two Raman 
beams. Example recoil frequencies En/h are given in 
table (|). 

We will assume throughout that rj < 1. This is the 
regime in which the ion trap processor is used in prac- 
tice, both because it facilitates the initial cooling to the 
ground state of motion, and because the processor then 
runs faster, as we will discuss. 



qubit basis 



transition 



recoil frequency (kHz) 



Beryllium 9 hyperfine ground states F=2 F=l 

Calcium 40 Zeeman ground states S±i/2 

Calcium 40 S1/2 ground and D 5 / 2 metastable state 



313 nm Raman 
397 nm Raman 
729 nm single-photon 



452 
63 
4.7 



TABLE I. Relevant parameters for the implementation of quantum information in trapped ions: wavelength, recoil frequen- 
cies for qubit candidate ions, assuming 45° angle between laser wave vector and the z axis. 



The Cirac-Zoller method adopts the internal state of 
each ion as a qubit (we restrict attention to two inter- 
nal states and thus one qubit per ion). Gates between 
qubits are obtained via excitation of the common vibra- 
tional mode, and measurement of the state of one or more 
qubits is by observing fluorescence. The method uses 
the fact that arbitrary single-qubit rotations, combined 
with any one 2-qubit gate such as controlled-not ( C X) 
or controlled-rotation ( Z), between arbitrary pairs of 



qubits, form a universal set |lCj[| . That is, any unitary 
transformation can be decomposed into elements from 
the set. In the trapped ion system, a single-qubit rota- 
tion is a transition in the internal state of a single ion, 
driven by a laser pulse. A C X between the internal states 
of any pair of ions A and B is achieved by a swap (S) 
operation between the internal state of A and the vibra- 
tion, followed by X between the vibration and the inter- 
nal state of ion B (with vibration as control, ion internal 



2 



state as target) followed by S again on ion A. Therefore 
the only operations that will concern us are single-qubit 
rotations, and S and C X between a single ion and the 
vibration. All of these are achieved by excitation of a 
chosen transition. The unwanted off-resonant excitation 
of other transitions are the main subject of this paper. 

The strength of the ion-laser interaction is 
parametrized by the Rabi frequency, given by 



Qnm = (n| exp (ir?(a t + a)) |m) f2f reo (1) 
= C nm Vl (2) 

where f2f rcc is the Rabi frequency for a free ion, 
the states |n) are vibrational energy eigenstates, Q = 
exp(— 77 2 /2)f2f rcc and the factor C nm is given by the fol- 
lowing equation (3). Values of C nm are listed in table 
(Q) for the low-lying vibrational levels. 



min ("V") f_1V' 2j 

' nm = V^!(^) |/ - m| E ^TFXl m I ( r-^T- 

~^ ]\3 + \ n ^ m\)\(pxm[m, n) — j)\ 



Cnm 





1 


2 


3 




1 

2 
3 


1 

irj 

-V 2 /V2 
-in 3 Mi 


ir] 

(l-» 2 ) 
iV2 V (l - 7)72) 

-y/3/2r, 2 (l- V 2 /3) 


-V 2 /V2 
iV2 V (l - 77 2 /2) 
(l-2n 2 +r7 4 /2) 
iV3r)(l - v 2 + V 4 /®) 


-irj 6 /V& 
-^/3/2 V 2 (l- V 2 /3) 
iV3r)(l -v 2 + " 4 /6) 
1 - 3t7 2 + 3ry 4 /2 - rf /6 



TABLE II. Matrix element for vibrational-state-changing transitions. 



The single-qubit rotations can be much faster than the 
two-qubit gates, because the energy separation Tiujq of the 
internal energy levels is much larger than that of the vi- 
brational levels, and by choosing a laser beam direction 
k z = (single photon transition) or k lz — k 2z (Raman 
transition), the Lamb-Dicke parameter can be made to 
vanish during single-qubit operations, which means these 
operations do not couple to the vibrational state (C nm 
becomes 8 nm ). Therefore, for An = 0, VL can be large 
compared to uj z without causing off-resonant excitation 
of An transitions. Therefore, the speed of the ion 
trap processor is limited by the S and °X gates. 

The S gate is achieved by a ir pulse on either the first 
red vibrational sideband, that is, at frequency u>q — lo z , or 
the first blue vibrational sideband at frequency ujq + oj z . 
The duration of the S gate is therefore 



(4) 



where the gate becomes exact in the limit Q — * 0. The 
choice of red or blue sideband can be dictated by exper- 
imental convenience. The two are equivalent in terms 
of their quantum computational effect, since the logical 
operations produced can be made identical simply by re- 
labelling the states (i.e. changing which internal state of 
the ion is called 0, and which is called 1). We will treat 
the red sideband throughout our theoretical discussion, 
but the results will apply equally to blue sideband excita- 
tion. In fact, we used a blue sideband in the experiments 
described in section III. 



We will consider two methods for the X gate. The 
first is that described in the original proposal of Cirac 
and Zoller, where X is obtained by single-bit rotations 



(Hadamard gates) combined with C Z, and C Z is obtained 
from a 2tt pulse on the first red sideband (or blue, de- 
pending on the relative positions of the levels) of an aux- 
iliary transition in the ion, i.e. laser frequency w aux — ui z , 
gate duration 



T - 2W 
TC1 ~ W 



(5) 



The other method is that of Monroe et al. Q, who 
proposed using a 2ni7r pulse on the carrier (frequency 
luq) where m is an integer and n 2 = l/(2m). This is es- 
pecially useful if an auxiliary transition is not available. 
When the vibrational state is |n = 0), this drives an in- 
teger number of Rabi oscillations of the internal state, 
while if the vibrational state is |n = 1), this drives a half- 
integer number of Rabi oscillations, since Cn = 1 — n 2 — 
(2rn— l)/(2m). The method requires 77 to be fixed in the 
experiment to one of the special values 1 / >/2m, which is 
easily done in practice by adjusting the trap confinement 
and/or laser beam direction. The duration of a C X gate 
by the Monroe method is 



T, 



C2 



2to7t 

n 



TV 



(6) 



Since the ion trap is operated with r\ < 1, it might be 
thought that Tq2 is necessarily greater than Tqi- In fact, 
this is not necessarily the case, since these gate times are 
limited by the maximum allowable Rabi frequency fi, and 
this can depend on the type of gate. 

Sorensen and Molmer |ll],|l2| have proposed a general 
method to implement gates such as controlled-not us- 
ing bichromatic laser fields, achieving good fidelity even 
when the vibrational degrees of freedom are not in their 
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ground state. The speed limitations of this method have 
recently been re-examined Jl3| , so we will not discuss it 
here, other than to say the method appears to be useful. 



B. Solution of time-dependent Schrodinger equation 

Our aim is to find the maximum switching speed of 
the processor. It is seen from equations @), and (|^) 
that this is determined by the maximum allowable Rabi 
frequency during the gates which operate on the vibra- 
tional state. The Rabi frequency cannot be arbitrarily 
large, because in the limit f2 3> to z , the evolution would 
be independent of the vibrational state. 

In order to find the maximum Rabi frequency, we need 
to solve the time-dependent Schrodinger equation for the 
system, without making the approximation of small Q. 
We assume a two-level ion. This means we will not ex- 
plicitly examine the type of C X gate which uses an auxil- 
iary level (equation (^|)), but in any case this operation is 
closely related to the S operation which we will examine, 
so results for the rate of S will apply to this type of C X 
apart from the factor 2. We will study S and the Monroe 
C X gate. 

The Hamiltonian for the ion that is being illuminated 
by the laser during a given gate is 



H = H + Hi 



(7) 



where Hq is diagonal, with diagonal elements given by 
H(0, w Zl 2lj z , ■ ■ ■ , wo, + u z , u>q + 2uj z , • • •), and the inter- 
action Hamiltonian 



H, 



9. 



Ce lULt 



(8) 



where C is the matrix having elements C nm . Note that 
the only approximation so far is to ignore the other vi- 
brational modes. 

We first adopt a frame rotating with the laser 
frequency: \fp(t)) = U\ t/j(t)) , where [/ is a 
diagonal unitary matrix, with diagonal elements 
exp(i{— 1, — 1, — 1, • • • , 1, 1, 1, • • -}uiLt/2). The states 

|V>(£)) satisfy the Schrodinger equation ih(d/dt) ?/>\ — 

H \ij>j, where H = ih(dU/dt)W + UHW is now a time- 
independent Hamiltonian. We can therefore write the 
solution to the Schrodinger equation 



\m)=e 



_ -im/h 



-iVHV'h 



/h V\ m) (9) 



where V is the matrix of eigenvectors of H , so that VHV^ 
is diagonal. 

The quantity of interest, from the point of view 
of quantum information processing in the trap, is the 
final state expressed as a superposition of computa- 
tional basis states. If at t = the computational 



basis states are |u(0)) = \u(0)), then at other times 
they are exp(—iHot/U) \u(0)) since then the coefficients 
(u(t) \i>(t)) do not evolve in the absence of gate opera- 
tions. Therefore we would like to calculate 



'(0)1 ■ 



-i(H-H Q )t/h 



(10) 



Define P = exp(—i(H — H^t/K), and let G be the 
precise unitary operator for the intended gate. Then the 
degree to which the laser pulse produces the intended 
gate is given by the overlap between the final state and 
one which would be obtained from G. To indicate the 
degree of imperfection of the laser pulse operation, we 
therefore calculate 



f min = mm\(iP\G^P \il>) | 



(11) 



and define the imprecision to be e = (1 — /mm) 1 ^ 2 - 

Let us comment on whether or not e ^ represents 
imperfection in the system. We are assuming no imper- 
fection in the sense of unknown dynamics (e.g. laser in- 
tensity and frequency noise) . Therefore as long as we can 
calculate numerically the effect of the laser pulse opera- 
tion, we have accurate knowledge of the expected state 
of the system: the fact that the laser pulse gate differs 
from any particular "ideal" gate is not a source of any 
imprecision at all. However, the intention is to use the 
ion trap as a quantum processor, to do quantum calcu- 
lations which we lack the computing power to simulate 
classically. Putting together many laser pulses, we are 
therefore assuming we cannot predict the effect on the 
whole quantum process of having e ^ in all the oper- 
ations. Therefore e must be regarded as imprecision in 
the device, which must be minimized. 

The question of driving complicated evolution in a pre- 
dictable manner is central to other techniques such as 
NMR spectroscopy, and there exist methods to combine 
pulses causing different types of driven rotation in order 
to undo the effect of certain terms in the Hamiltonian. 
Equivalent methods can almost certainly be found for 
the ion trap system, but in any case operations close to 
ideal ones will remain the best starting point for any such 
method. 



C. Performance without correction for light shift 

In this section, we will examine the operation of the 
processor in the regime 77 < 1, in the case that the light 
shifts (a.c. Stark shifts) are not taken into account when 
choosing the laser pulse frequency, phase, and duration. 

First consider the S gate: a 7r pulse on the first red 
sideband. We set = ujq — lo z in (II B), and exam- 
ine the eigenvalues of H . The separations between the 



eigenvalues enable one to deduce the Rabi flopping fre- 
quencies and the light shifts in the system. The Rabi 
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flopping frequency on the |n = 0) «-> \n = 1) transition 
is r]Q, as expected, which confirms that the gate time is 
Ts = Tr/(r]tt) as in equation (|]). 

The frequency of the transition is shifted by 



A W = i|l 



>£ 

2 



1 



(12) 



This is primarily the light shift caused by the presence 
of An = transitions, which are off-resonant by uj z . 

As was noted by Wineland et al. fl4f| , the primary er- 
ror in the gate operation is caused by the light shift Auj, 
which can be significant compared to the Rabi flopping 
frequency 77!). To understand its effect, we now model 
the transition of interest, that is \g, n = 1) <-> |e, n = 0), 
as a two-level system, driven by a 7r pulse off-resonant by 
Auj. The only term in the propagator of first order in Auj 
is i{Auj/a) sin(at/2) where a 2 = rj 2 Q 2 + Auj 2 . Putting 
rjQt = 7r, we obtain 



iAuj 



n 



2t]uj z 



(13) 



Therefore, to attain a given gate precision e the Rabi 
frequency must satisfy fl < 2tr\uj z , which gives 



1 



r?H 2e 2 
— < — 77 UJ Z 

7T 7T 



= 4e 



E 



R 

Nh' 



(14) 
(15) 



The result thus has the simple form that the processor 
speed for swap operations is limited by the iV-ion recoil 
frequency. 

Next, consider the ^Xgate using the special Lamb- 
Dicke parameter method of Monroe et al., that is, a 
2m7r pulse on the carrier with 2m = 1/n 2 . Substitut- 
ing ujl = wo into (II B), and examining the eigenvalues 
of H, we find the transition frequencies are shifted by 



Auj 



uJ z 



(16) 



Since Auj <C fi, the light shifts for this gate are much less 
significant. Modelling each driven transition as a two- 
level system, using the same method described above for 
the S gate, we find that the effect of Auj is e ~ r] 2 Q,/uj z . 
This is small compared to the effect of off-resonant ex- 
citation of An = ±1 transitions. These are driven in 
the regime where the detuning is large compared to their 
Rabi frequency, leading to a propagator with terms of 
order ^J2rfljuj z , so we find 



Setting < euj z /\/2r], we find 



(17) 



1 ^ F> \ Er Uz 



n<t:uj z 



(18) 



It might appear from this that the gate speed can in- 
crease without limit by using a tight trap, but this is a 
false impression because the expression is only valid in 
the regime f2 <^ uj z , which means it is only valid for a 
gate rate small compared to i] 2 uj z /n ~ Er/{NK). 



D. Gates with correction for light shift: Red 
sideband 



The light shift is readily corrected for by tuning the 
laser frequency to accurate resonance with the light- 
shifted transition frequency We will now describe 
the correct laser pulse frequency, phase, and duration, 
and derive an expression for the imprecision caused by 
the primary remaining unwanted effect, which is off- 
resonant excitation. We find that the recoil frequency 
still strongly influences the processor speed. 

We have studied this problem by numerically evalu- 
ating expressions (|l^) and ([ll]). Only vibrational levels 
with n < 4 were included, thus yielding an 8 x 8 propaga- 
tor matrix. When further vibrational levels were included 
in the calculation, the results were not significantly af- 
fected. 

Wc find that the best implementation of a swap gate is 
obtained when the laser frequency is tuned to the light- 
shifted frequency of the first red sideband: ujl = uJo — 
uj z + Auj where Auj is given by equation (|l^). The gate 
time is still as in equation (^). With no further adjust- 
ments, this produces a propagator close to U^S where S 
is a perfect swap operator, and is diagonal with diag- 
onal elements exp(z{l, 1, 1, • • • , — 1, —1, —1, • ■ -}AluTs/2). 
Therefore, to produce the intended gate, the system must 
be further evolved by an application of U\. This is equiv- 
alent to applying subsequent laser pulses with a corre- 
spondingly adjusted phase. Therefore the swap gate is 
completed by changing the phase of the oscillator used 
in the experimental apparatus to keep the laser in step 
with the atom's internal resonance (at frequency ujq), by 
an amount Aip — AujTs- 

This adjustment permits the swap gate rate to exceed 
the recoil frequency. The term in the system propaga- 
tor P which now gives the dominant contribution to the 
imprecision e is off-resonant excitation of An = transi- 
tions. We expect the amplitude for this process to be of 
:. Our numerical calculations confirm 



order n/^/W 
this, giving 



^/2(n 2 +uj 2 ) V2uj z 



(19) 



for the remaining imprecision in the adjusted swap oper- 
ation. Therefore the gate rate is now limited by 
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which is 2^/2 e times the geometric mean of the iV-ion re- 
coil frequency and the vibrational frequency of the chosen 
normal mode. 

If an auxilliary transition is available, the above argu- 
ment will hold for the C Z gate implemented by a 2ir pulse 
to the first vibrational sideband of the auxilliary transi- 
tion, the only difference is that the gate rate at given e 
is half that for the swap gate (T C i = 2T S ). 



recoil frequency, independent of the trap tightness, but 
the gate becomes more precise as the trap gets tighter 
(e ~ i]/V2, expression (|l7|)). 

One way to use the carrier excitation method to 
achieve faster X gates, is to use higher vibrational levels 
such as \n= 2) (in conjunction with |n = 0)). However, 
this would mean that the S gate would have to drive the 
second red sideband, making it slower (the expected lim- 
itation at small r\ being off-resonant excitation of An = 
transitions). Therefore only small gains, if any, are avail- 
able by this route. 



E. Carrier 



F. Sensitivity to laser intensity fluctuations 



The Monroe method X gate is executed by a laser 
pulse at the carrier frequency. The speed limit given by 
equation (|l^) cannot be exceeded, since it is caused by 
off-resonant excitation that cannot be avoided. However, 
we are interested in the behaviour outside the region of 
validity of (|l8|), that is, when ~ u> z . 

When 7] is small but fl is not, the precision of the 
operation need not be limited by off-resonant excitation 
of An = ±1 transitions, since their rate is low. The 
source of imprecision instead comes from the fact that 
the dynamics of the system can no longer be pictured 
as separate Rabi flopping on the n = <-> transi- 
tion and the n = 1 1 transition. Rather, the inter- 
action Hamiltonian drives the set of vibrational levels 
n = 0, 1,2, • • • ~ Vt/ui z as a single entity. The result is 
that the evolution from an initial state \n = 0) is almost 
the same as the evolution from an initial state \n = 1). 
Therefore the Rabi flopping must be driven for a longer 
time before a C X gate is obtained, and the gate speed 
ceases to increase with f2. The rate obtained at fl ~ lo z 
was found to be the maximum. Recalling the expression 
(||) for Tc2 , it is seen that this maximum is approximately 
the recoil frequency. 

In the region f2 ~ ld z , the gate can be optimized sim- 
ply by adjusting the pulse duration. The correct pulse 
length is not precisely ir/(r] 2 Q), but differs from this by 
an amount of order 7r/(ni7). We have not found any sim- 
ple expression for this adjustment, we surmise that this 
is because it is a complicated function of all the light 
shifts in the multi-level system. To run a processor in 
practice, the correction to the gate time can either be 
calculated numerically, or measured experimentally By 
this adjustment, we were able to extend the region of 
validity of expression ( [l8| ) up to f2 ~ 0.5lu z . Operating 
at this limit, the gate rate is always equal to the TV-ion 



Although our main purpose is to study limitations im- 
posed by the unavoidable properties of the system, we 
comment here on the sensitivity to one source of techni- 
cal noise, namely laser intensity fluctuations, since these 
partially limit our experiments described in section III. 
The laser drives Rabi flopping at frequency ryJ7 for the 
Cirac-Zoller swap gate, and at frequency fl for the Mon- 
roe gate. If, owing to intensity fluctuations, f2 is impre- 
cise by A17 then the action of either gate will be impre- 
cise by e ~ TsnAVl/ir and e ~ Tc2Aft/ir respectively. 
We therefore find AQ/fi ~ e for the swap gate (also for 
the closely related Z gate), whereas AQ/Cl ~ n 2 e for 
the Monroe gate. 



G. Allowance for all the normal modes 

The only approximation we made was to neglect all 
but one of the normal modes of oscillation of the ion 
string. When we relax this assumption, the main fea- 
tures of the more general problem can be understood in 
terms of the Rabi frequencies and the normal mode fre- 
quencies. To keep the notation clear, we will describe 
a case of three normal modes, which we will label x, y 
and z, but the discussion readily generalises to more nor- 
mal modes. In our experiments described in section III 
there is a just a single ion, so the modes of interest are 
indeed vibration along the principle axes of the trap, but 
the modes here labelled x, y and z could equally refer to 
different normal modes of the z-oscillation of a three-ion 
string, in the case where the two other directions of oscil- 
lation are frozen out owing to the linear geometry (that 
is, they have much higher frequency, and are prepared in 
the ground state). Equation (0) is replaced by 



n(n x ,m x )(n y ,m y )(n z ,m z ) 



(21) 



G 



in an obvious notation, where we now redefine il = 
exp(— (r) x + T)y + 77^) /2)S7f ree - There are two regimes to 
consider. In the case that all the modes are prepared 
very close to the ground state, the main effect is that the 
light shift is no longer that given by equation (|l2|), but 
one given by a sum of light shifts related to all the possi- 
ble transitions. This is easily compensated by a suitable 
adjustment to the laser frequency. The limit on precision 
therefore remains that due to off-resonant transitions. 
For excitation on the red sideband (n z = 0, m z = 1), 
in addition to carrier transitions excited off-resonant by 
u> z with Rabi frequency tt, there are also sideband tran- 
sitions excited off- resonant by oj x — u> z and uj y — lo z , with 
Rabi frequencies r) x £l and rj y fl respectively. Since the 
latter contribute at a higher order in r\ compared to the 
carrier term, they are negligible unless two modes are 
close in frequency. For a linear ion string the separation 
of the lowest mode frequencies is large enough that the 
additional terms make little contribution. Therefore our 



previous discussion, including equations (|19|) and (|2C 
remains valid. In the experiments described in section 
[II, on the other hand, to y — lo z ~ 0.025 lj z so unwanted 
excitation of the y sideband has to be taken into account. 

For the carrier excitation, in addition to the z side- 
bands excited off-resonant by u> z with Rabi frequency 
r) z fl, there are also sideband transitions excited off- 
resonant by u> x and to y with Rabi frequencies (1 — rj^)r) x Q 
and (1 — rji)riyVL respectively. These contribute at the 
same order in 77 as the z sidebands, and result in a change 
to equation ( |l8| ) by a numerical factor of order unity. 

Our experiments were carried out with only one mode 
cooled to the ground state. The other two modes 
had roughly thermal distributions with mean vibrational 
quantum number of order 10. The result of this is to 
blur all the Rabi frequencies and light shifts, since each 
time a given z transition is excited, the Rabi frequency 
depends on the x and y vibrational quantum numbers: 



ft(n x ,n x )(n y ,n y )(n z ,m z ) ~ (l - n x rj£) (l - n y r) y ) C, 



n. 



(22) 



We are able to fit our experimental results by calcu- 
lating the evolution for each value of n x ,n y , and then 
averaging over a thermal population of the x and y vi- 
brations. 



H. Summary 

In conclusion, we find that the swap gate and the 
Cirac-Zoller type of controlled-phase gate can be made 
faster, at a given level of precision, by making the trap 
tighter, but the speed increase is only in proportion to 
u>z . This is in agreement with earlier work but here 
we have added a precise quantitative statement of the 
trade-off between speed and precision. Some previous 
studies have adopted T$ ~ tt/w z for the purpose of mak- 
ing rough estimates, but we find the error thus introduced 
is significant; for example, it leads to an overestimate of 
the switching rate by two orders of magnitude in p5[ . 
By exciting the ion in the node of laser standing wave, it 
is possible to avoid (or greatly reduce) the off-resonant 
carrier excitation , therefore T$ ~ 7t/oj z may be avail- 
able, but this involves a significant extra complication of 
the experiment which may render it unfeasible in practice 

The "magic Lamb-Dicke parameter" method of Mon- 
roe et al can not be speeded up indefinitely, it has a 
natural speed limit given by the recoil frequency. 

There are advantages in having the ions spaced by a 
micron or more, in order to allow laser addressing of one 
ion at a time, and/or resolving the fluorescence from dif- 
ferent ions with appropriate imaging optics Jl7[ ). If we 
set a limit s for the smallest permissible seperation of 



the closest ions in the string during processing, we ob- 
tain Ail,!. 



W, '°» < 47re Afs 3 iV 1 - 71 



(23) 



Putting this in equation (ppf), and adopting the breath- 
ing mode u> z = \/3wz,cm, we obtain for the fastest swap 
gate rate, at given precision and ion spacing, 



1 



T S ~ ' WW I h J MW/4JV0.93- 



(24) 



If we approximate the N dependence as TV -1 , and set 
s = 10 A where A is the laser wavelength, then we arrive 
at a "gate time per ion" which depends only on the choice 
of ion and transition, at given e. This quantity, Tg/N, 



is given in table [II for some example transitions, for the 
case of 99% gate fidelity (i.e. e = 0.1). 



ion 


mass 


wavelength En/h (kHz) 


gate time per ion 


Beryllium 


9 


313 nm 452 


1.26 /is 


Calcium 


40 


397 nm 63 


5.6 ^is 


Calcium 


40 


729 nm 4.7 


34 /is 



TABLE III. Minimum gate time per ion for 99% fidelity 
SWAP gates, for three example transitions. The gate time per 
ion, Ts/N, is given by equation (pif), for the case e = 0.1, 
s = 10A ~ 3,4, 7 nm. 
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III. EXPERIMENTS 



In this section, we describe our experimental investiga- 
tions, carried out in conditions where the Rabi frequency 
^sideband = ^carrier of the sideband excitation is signif- 
icantly above the recoil frequency. We thus are able to 
investigate the regime where the speed of execution of 
the gates determines the fidelity of operation. 

We use a single trapped 40 Ca + ion to demonstrate the 
principle of quantum information processing. The elec- 
tronic \Si/2, TTi = —1/2) ground state and the metastable 
l-^5/2j m = ~ 5/2) state (Is lifetime) are used to imple- 
ment one qubit. We apply a 4 Gauss bias magnetic field 
to lift the degeneracy of sublevels in the ground and ex- 
cited state manifolds. The qubit can be coherently ma- 
nipulated by laser light at 729 nm (see table |). Two ba- 
sic operations are demonstrated. First there is the single 
qubit rotation which only affects the individual internal 
electronic state, and secondly we perform the swap oper- 
ation which entangles the electronic state of the ion and 
its vibrational state. We did not investigate the Monroe 
et al. X gate because our experiments involved small 
Lamb Dicke parameters of order 0.045. In this regime the 
Monroe method is ~ 500 times more sensitive than the 
Cirac-Zoller method to laser intensity noise and thermal 
populations in spectator modes, and we found it to be 
unworkable in these experiments. 

We store the ion in a conventional spherical- 
quadrupole Paul trap with trap frequencies uj x .y,z = 
27r(4.0, 1.925, 1.850) MHz. To acquire each experimental 
data point, at given values of the parameters, we run a 
series of 100 cycles. A single experimental cycle is made 
up of 5 stages, which are (i) Dopplcr cooling, (ii) side- 
band cooling, (iii) coherent driven evolution, (iv) obser- 
vation of fluorescence, (v) deshelving. We then record 
the fraction Po of the 100 cycles in which fluorescence 
was observed in the penultimate stage. 

In more detail, the 5 stages of a cycle are as follows. 
Dopplcr cooling is performed on the \S\12) to \P\/2) 
397 nm electric dipole transition. The electronic ground 
state is then prepared in a pure \S\/%,m — —1/2) state 
by optical pumping. The ion is cooled to the vibrational 
ground state of the oj z mode by applying sideband cool- 



ing on the \S\12 



-1/2) 



\D. 



5/2i 



-5/2) elec- 



tric quadrupole transition at 729 nm. A full description 
of the trap and the cooling procedure is given in Ref. 
H . The wavevector k of the 729 nm radiation is nearly 
perpendicular to the y direction, and has an angle of 40° 
and 50° to the z and x directions. The corresponding 
values of r^ 9 " m are (0.04, 0.01, 0.045). In the qubit op- 
eration step, the \S\/2,m = —1/2) «-> \D 5 / 2 ,m = —5/2) 
transition is excited with a laser pulse of well controlled 
frequency, intensity, and timing. Then, by monitoring 
fluorescence at 397 nm under laser excitation, we detect 
whether a transition to the non- fluorescing state D 5 / 2 



occurred. The scheme allows one to discriminate be- 
tween the internal states of the ion with an efficiency 
close to 100% JT^-pOt- In our experiment, the discrim- 
ination effiency is approximately 99.8%, limited by the 
lifetime of the metastable state pl[ . 

Finally, the ion is again repumped (deshelved) from 
the D5/2 level to the electronic ground state via the P 3 / 2 
level. The fraction Po of cycles in which fluorescence is 
observed indicates the population of the D 5 / 2 level after 
the coherent qubit operation step. 



A. Carrier excitation 



As discussed in section [I A , the single qubit rotation 
can be driven fast and without loss in contrast, even in 
the range where O exceeds the trap frequencies u> : 



x,y,z ) 



as 



long as the Lamb-Dicke parameters are sufficiently small. 
To reach this limit in the case of a single-photon tran- 
sition under travelling wave excitation, the wavevector 
k of the exciting light field should be perpendicular to 
the chosen vibration direction, and to the direction hav- 
ing next smallest vibrational frequency. The remaining 
vibrational frequency must be large compared to CI. In 
our case the laser is nearly perpendicular to one of the 
weak confinement directions, but not the other, so there 
is some unwanted excitation of the motion. A limit to ex- 
ecution speed is set by the available laser power to drive 
the qubit transition. However, Rabi frequencies as high 
as a few MHz can be achieved even in the case of dipole 
forbidden quadrupole transitions. 




40 60 80 100 

Time (|is) 

FIG. 1. a) Rabi oscillations on the carrier. We excite the 
An = transition with a laser power of 100 mW and observe 
a SI = (2n) 1090 kHz Rabi frequency. A contrast of better 
than 95% is reached for n and 2tt pulses, b) Rabi oscillations 
on the blue sideband at 4mW excitation power. At this low 
power the contrast for ir and 2tt pulses is 95%, respectively 
80%. 



8 



Rabi oscillations on the An = carrier transition at 
a frequency ^carrier = 27r 1090 kHz are shown in Fig. 
la. For this, lOOmW of light are focussed into a waist 
of 30/xm at the position of the ion. We observe a con- 
trast of more than 95%, the contrast being mainly lim- 
ited by the thermally distributed axial vibration mode 
(u) x = 2tt 4.0MHz), which causes a spread in Rabi fre- 
quency owing to the n-dependence as indicated in equa- 
tion (|2^). We find that our results are consistent with 
an average over a thermal phonon distribution in the x- 
oscillator having (n) = 12(2). 



l D 5/2' 0> l°5/2. 1> 




|S 1/2 ,0> |S 1/2 , 1: 



FIG. 2. Relevant energy 
levels: \Si/ 2 ,n = 0), 
\D S /2,n = 0), \S 1 / 2 ,n = 1), 
and 

\Ds/2,n — 1). The solid ar- 
row indicates the blue side- 
band excitation, the dotted 
one the off-resonant carrier 
excitation. 



B. Sideband excitation 



The second building block is the sideband excitation, 
involving a An ^ transition, here performed on the 
z vibration. For convenience, we adopted the first blue 
rather than red sideband. The equivalence of blue and 
red sidebands for these studies was mentioned in sec- 
The laser frequency was red-detuned from the 



II A 



tion 

blue sideband resonance (measured at low laser power) to 
compensate for the a.c. Stark effect. The accuracy of the 
An = +1 operation is therefore limited by off-resonant 
excitations on the An = transition, as discussed in sec- 
tion II D , see also Fig. 2. We indeed observe off-resonant 
carrier excitation, visible as fast oscillations with a fre- 
quency near uj z on top of the £l S ideband — 2tt 48 kHz Rabi 
oscillation, see Fig. 3. Note that this frequency is far be- 
yond the recoil frequency, and thus the contrast shrinks 
to 75% for a n pulse. 
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Time ( |is) 

FIG. 3. Rabi oscillation on the blue sideband at lOOmW excitation laser power. See text for details. 



If the blue sideband lu z is excited at a Rabi frequency of 
^sideband = 27r7.4kHz (Fig. lb), the off-resonant carrier 
excitation is no longer detectable. A contrast of 95% for 
a 7T pulse and 85% for a 2n pulse is measured, while equa- 
tion (|l9| ) predicts e = 0.063 and hence a fidelity 99.6% 
for the 7T pulse. Under low-power excitation, we do not 
observe any light shift. 

The experimental data can be fitted by numerically 
solving the Schrodinger equation in the truncated ba- 
sis of the 4-level system consisting of \Si/ 2 ,n = 0), 
I-D5/2! n — 0), |<Si/2j n = 1), and \D 5 / 2 , n = 1) states (see 
Fig. 2) , and averaging over a thermal distribution of pop- 
ulation in the other vibrational modes. Three parame- 
ters in the calculation were independently measured: the 
Rabi frequency on the carrier £1 = 27iT090 kHz, the Lamb 
Dicke parameter rj z , and the trap frequency lu z — 2tt 1850 



kHz. The detuning of the laser field 8 is varied for 
optimum contrast and we find excellent agreement be- 
tween the data and a numerical simulation for a value of 
<Sth/27r= -375 kHz (Fig. 3). The discrepancy of 125 kHz 
between 5 t h and the experimentally determined value 



>/2tt 



-250 kHz is probably caused by light shifts 



due to the D 5/2 - P 3 / 2 , S 1/2 - P1/2, and S 1/2 -P 3/2 tran- 
sitions. The sum of the calculated light shifts from the 
dipole transitions is approximately 70 kHz. 

A third contribution to the overall light shift is ex- 
pected from the other radial (y) mode at 1.925 MHz 
(5oj/2it=75 kHz). Each level of the y- vibrational ladder is 
light-shifted, resulting in a change in the (n y , n y )(l z , Z ) 
z-sideband resonance frequency by approximately (r\ v ■ 

^carrier 

) 2 (n y + l)/(2(uy - uo z )) = {n y + 1)(2tt) 0.8kHz. 
Assuming a mean vibrational quantum number of 25, a 
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20 kHz spread in shifts is expected. To model this sit- 
uation we solved the Schrodinger equation for different 
y-oscillator occupations and averaged the results over a 
thermal distribution (see Fig. 4) . The loss of contrast at 
fit = 7r is dominated by this thermal averaging, rather 
than laser intensity and magnetic field fluctuations. 

To compare the experimental findings with an opti- 
mum 7r pulse, for a given pulse length, we investigated 
numerically the parameter space and found that an op- 
timum 7r pulse (92% contrast of Rabi |S f _ 1 / 2 ) l-D-5/2) 
oscillation) could be achieved for a detuning of Sth/2n = 
—355 kHz. We conclude that the detuning chosen in the 
experiment was off by 20 kHz. Note, that the electron 
shelving detection method only measures the internal 
state occupation, not the vibrational one. Hence, the 
measured contrast does not directly yield the fidelity. 
Equation (|l9| ) predicts e = 0.41 for this case, and hence a 
fidelity of 83%. The 20 kHz deviation from the ideal de- 
tuning leads to a 5% loss in contrast and the thermally 
distributed y-vibration mode with (n) = (25) further re- 
duces the contrast to 75%. The corresponding fidelity 
was near 64%. For the details of the ir pulse dynamics, 
see Fig. 4. 



ments of the precision of the operations, taking into ac- 
count the complete Hamiltonian including all the vibra- 
tional states and the off-resonant coupling terms. These 
statements are given in equations ([l7]) , ( |l8| ) for the Mon- 
roe et al "magic Lamb Dicke parameter" gate, and in 
equations ([TJ), @) for the Cirac-Zoller swap gate. If 
we set a limit of a fixed number of laser wavelengths for 
the spacing of the ions in a linear trap, then we arrive at 
a gate time per ion, for the swap gate, which depends 
only on the precision and the choice of ion and transition. 
This time is given in equation (E3) and table III. 



In our experimental studies we have demonstrated a 
contrast of 75% for a 1/Ts = 100 kHz swap operation, 
and 95% for a I/T5 = 14 kHz swap operation on a single 
trapped ion. These gate rates are respectively 21 and 3 
times the relevant recoil frequency. In both cases the con- 
trast could be significantly improved by cooling a further 
vibrational mode. 




Time (pis) 

FIG. 4. n pulse on the blue sideband, detailed view. The 
data of Fig. 3 are plotted together with a fit taking into ac- 
count a thermally distributed radial y-mode (see text for the 
set of parameters) . The dotted line gives a prediction for op- 
timized parameters, where all relevant vibration modes are 
assumed to be cooled to the ground state, and the detuning 
is set for maximum contrast. A contrast of 92% is predicted 
under these optimum conditions, and the fidelity 1 — e 2 given 
by equation ( |ig| ) would then be 83%. 

To conclude, we have presented a theoretical discussion 
of two types of quantum gate which couple the internal 
ion state with its motion, and an experimental study of 
one of these. Our calculations lead to quantitative state- 
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